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ABSTRACT Apparent competition between species 
is believed to be one of the principle driving forces that 
structure ecological communities, although the precise 
mechanisms have yet to be characterized. Here we de- 
velop a model system that isolates phage-mediated inter- 
actions by neutralizing resource competition using two 
genetically identical B. bronchiseptica strains that dif- 
fer only in that one is the carrier of a phage and the 
other is susceptible to the phage. We observe and quan- 
tify the competitive advantage of the bacterial strain 
bearing the prophage in both invading and in resist- 
ing invasion by bacteria susceptible to the phage, and 
use our measurements to develop a mathematical model 
of phage-mediated competition. The model predicts, 
and experimental evidence confirms, that the competi- 
tive advantage conferred by the phage depends only on 
the relative phage pathology and is independent of other 
phage and host parameters. This work combines exper- 
imental and mathematical approaches to the study of 
phage-driven competition, and provides an experimen- 
tally tested framework for evaluation of the effects of 
pathogens/parasites on interspecific competition. 

INTRODUCTION Pathogen-induced damage to 
hosts, commonly observed as infectious disease, has 
been extensively investigated in humans and other an- 
imals 0,12. A patho gen can also confer a competitive 
advantage to one of two competing species through the 
process known as apparent competition [1, ^ |E l^i IE 111 ■ 
In both scenarios pathogens appear to drive the evolu- 
tion of their hosts, exerting a selection pressure toward 
greater resistance ^^^tAcM- Increased resistance mech- 
anisms of the hosts, including those as complex as the 
adaptive immune response of higher eukaryotes, does not 
seem to confer freedom from infection, but offers sub- 
stantial advantage against other hosts more susceptible 
to the pathogen jOj- Pathogens are also driven to maxi- 
mize their fitness (a function of virulence and transmis- 
sion) although this is often a consequence of balancing 
virulence with transmission [Tol . Pathogen-mediated 
competition is an outcome of this ever-escalating arms 
race between the co-evolving hosts and pathogens. 

There are a few excellent illustrative examples from 
the laboratory and field studies of ecological assem- 
blage 0,0, |a I3l and even from the history of human 
diseases [i3,ll3|- However, due to the complexities origi- 
nating from dynamical interactions among multiple-hosts 
and multiple-pathogens, it is not always easy to single 
out and quantitatively measure the effect of pathogen- 
mediated competition in nature. In a system of bacte- 
ria and bacteriophage it is relatively easy to manipulate 



both host resistance mechanisms and pathogen virulence 
and thus this is one of the most suitable systems for the 
exploration of pathogen- mediated competition. In fact 
in recent studies 0, |3| laboratory communities of bac- 
teria and lytic bacteriophage have been used as model 
systems for phage-mediated competition between phage- 
sensitive and phage-resistant bacteria. However, in these 
studies resource- and phage-mediated competitions were 
strongly intertwined and the role of one of the most im- 
portant players in phage-mediated competition, lysogens 
(carriers of the phage), was not investigated. 

To address these concerns in an examination of 
pathogen-mediated competition, we established an in- 
fection system using bacteriophage, neutralized resource 
competition by having two genetically identical strains 
and using large nutrient excess, and then examined com- 
petition between bacterial strains that differ only in their 
sensitivity to this phage. Particularly we used Bordetella 
bronchiseptica, a causative bacterium of mammalian res- 
piratory disease, and its natural virus (BPP-1), a tem- 
perate phage that can either incorporate its DNA into 
the genome of B. bronchiseptica (lysogeny) or replicate 
itself and lyse the host bacterium (lysis). Here we demon- 
strate, both experimentally and theoretically, the exis- 
tence of a competitive advantage conferred by this virus 
in a bacterial population, using a system in which this ef- 
fect can be measured and quantified. Our results suggest 
that the lysogens are not only the source of phage dur- 
ing an infection process but can lead to fundamentally 
different dynamics in phage-mediated competition. We 
observe that the bacterial strain bearing the phage has 
an advantage in both invading and resisting invasion by 
bacteria susceptible to the phage and that the differential 
pathology on the two hosts is the sole variable that deter- 
mines the quantitative value of phage-mediated compet- 
itive advantage. The theoretical representation of these 
interactions should have broad application to pathogen- 
mediated competition at many levels. 

RESULTS The quantitative assessment of pathogen- 
conferred competitive advantage is usually hampered by 
the existence of direct competition over resources 3, 4, ^ 
. To address this concern we use as our host populations 
B. bronchiseptica strains that are genetically identical ex- 
cept for defined genetic changes. The wild type parental 
B. bronchiseptica strain RB50 (Bb) was used to gener- 
ate a strain (BbGm) that carries a gentamycin resistance 
marker shown not to affect expression of nearby genes 
and another (Bb::(/)) that is the carrier of, and is therefore 
resistant to, the temperate phage BPP-1 {(j)) 0,10]. To 
examine possible nutrient-dependent competition among 
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FIG. 1: The in vitro growth curves of B.bronchiseptica 
strains, (a) The three B.bronchiseptica strains, Bb::^ (open 
squares), BbGm (filled circles) and Bb (filled diamonds) have 
identical growth curves when grown separately. The straight 
line corresponds to a doubling time of 77 minutes, (b) Bb 
and BbGm grow without competition when co-cultured. 



these bacterial strains, we observed their in vitro growth 
rates. All three strains grew in nutrient rich medium 
with identical doubling time (77 minutes) over a range 
of cell densities from less than 1000 CFU/ml to over 10^ 
CFU/ml (see Fig. Bb and BbGm grew at the same 
rate when co-cultured, indicating that there is no effect 
of direct resource competition on their growth rates. The 
proportion of the final density of each strain was identical 
to the proportion at the start of the co-culture, indicat- 
ing neutral competition between these two strains under 
these conditions. 

We then examined phage-mediated competition using 
these B.bronchiseptica strains and the temperate bacte- 
riophage BPP-1. The interactions involved in this sys- 
tem are schematically represented in Fig.|2 All bacterial 
strains divide with a constant rate a and bacterial popu- 
lations grow with a density-dependent rate r. Susceptible 
bacteria (BbGm) become infected with a rate k, defined 
as the number of contacts between a phage particle and 
a host bacterium per unit time multiplied by the prob- 
ability of the host being infected upon contact.^ Upon 
infection the phage can take one of two pathways 18]. In 
a fraction P of infected BbGm, the phage replicate and 
then lyse the host after an incubation period 1/A, during 
which the bacteria do not divide Ts]. Alternatively the 
phage lysogenize a fraction 1 — P of their hosts, incorpo- 
rating their genome into that of the host. Thus the pa- 
rameter P characterizes the pathogenicity of the phage, 
incorporating multiple aspects of phage-host interactions 
resulting in damage to host fitness. The lysogens (Bb::0 
and BbGm::0) carrying the prophage grow, replicating 
prophage as a part of the host chromosome, and are (/>- 
resistant. Even though these lysogens are very stable [T3 
without external perturbations, spontaneous induction 



FIG. 2: Diagrammatic representation of phage- mediated 
competition between two bacterial strains. The phage {(p) 
is represented by a hexagon carrying a small thick line (0 
DNA). Bacteria are represented by a rectangle containing an 
inner circle (bacterial DNA) while the bacterial strains bear- 
ing the prophage (Bb::0 and BbGm::(/)) are represented by 
rectangles containing (j) DNA integrated into bacterial DNA. 
Gm stands for gentamycin marker. All bacterial populations 
(Bhwcj), BbGm and BbGm::(^) grow with identical density- 
dependent growth rate r. Lysogens are spontaneously induced 
with a rate 5. The rate of a single bacterium being infected by 
a phage particle is given by the infection-causing contact rate 
K. A is the infection-induced lysis rate, and P the probability 
of the phage taking a lytic cycle. Resource is accessed and 
consumed by all bacteria. 



can occur at a low rate 6, consequently replicating the 
phage and lysing the host bacteria. In general, both the 
number of phage produced (burst size) and the phage 
pathology P depend on the culture conditions J^] . This 
model predicts that when co-cultured, these two strains 
should compete through the phage, giving advantage to 
the lysogens. 

To experimentally observe phage-mediated competi- 
tion between two strains, Wo::(f> and BbGm were co- 
cultured in vitro. Since spontaneous induction of the 
phage from Bb::^ is inevitable but would occur at vari- 
able time points (see Appendix B), we added a small 
number (1,000 PFU/ml) of exogenous phage to syn- 
chronously initiate phage-mediated competition. 1,000 
CFU/ml of Bh::(t) and 1,000 PFU/ml of phage were 
added to a culture containing 10,000 CFU/ml of BbGm. 
As shown on Fig. Ol^a), the total BbGm concentration 
(both susceptible BbGm and BbGm::^) increased with 
the same initial growth rate as Bh::cf), but suddenly fell 
at about 8-12 hours postinfection. Within 24 hours the 
initial 1:10 ratio of Bb::0 to BbGm was reversed to ap- 
proximately 10:1, indicating a 100 fold relative increase 
in the proportion of Bb: :<?!). 

To observe the advantage of the lysogens in resisting 
invasion by the strain susceptible to the phage, 1,000 
CFU/ml of BbGm was added to a culture containing 
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FIG. 3: In vitro experiments (symbols) and numerical sim- 
ulations (lines) of (a) the invasion of the strain (Bb::(^) ex- 
ogenously and endogenously carrying the BPP-1 phage to 
the susceptible strain (BbGm), and (b) the protection of 
Bb::(^ against the invading BbGm. Symbols and lines rep- 
resent Bb::(/!> (open squares, thick solid line), BbGm::(/!> (open 
circles, thin solid line), and the total BbGm (filled circles, 
long-dashed line), respectively. The parameters used for 
the numerical simulations are x ~ 50, a = 5 /a — 0.1, 
P = X/a = 0.15, P = 0.98. 7 = SB{0)K/a = 0.01 for (a) 
and 7 = 0.002 for (b) where Sb{0) is the initial concentration 
of BbGm. 



10,000 CFU/ml of Bb::0 and 10,000 PFU/ml of phage 
in Fig.Otb). The initial 10:1 ratio of Bb::(/) to the total 
BbGm was amplified to approximately 1,000:1, indicat- 
ing again a 100 fold increase in the proportion of Bb::0 
to the total BbGm. 

We constructed a theoretical model of phage-mediated 
competition based on susceptible-infectious (SI) models 
describing the impact of directly transmitted pathogens 
on the dynamics of multiple hosts jl9jM,-21, 22, 23, 24]. 
We consider the dynamically interacting system of five 
subpopulations: (1) bacterial strain A (Bb::(/)) which 
bears the prophage and is thus resistant to the superinfec- 
tion by the phage, bacterial strain B (BbGm) which can 
be in one of (2) susceptible, (3) latent, or (4) lysogenic 
states, and lastly (5) the bacteriophage ((/>). Here a la- 
tent bacterium is one which is currently infected with the 
phage and will be lysed after the incubation period of the 
phage. We assume homogeneous mixing among all sub- 
populations, justified by vigorous stirring of the culture 
growth tube on an agitator. The time-evolution of each 
subpopulation is determined by the rates of incoming and 
outgoing flow in subpopulations, quantified by the rates 
given in Fig. [3 The five dimensional system of ordi- 
nary differential equations describing the bacteriophage- 
mediated competition is given by 
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where I{t), L{t) and S{t) are the concentrations of the 
infected, latent and susceptible bacteria at time t, N{t) = 
lA{t) + iBit) + Seit) + LBit), r{t) = a(l - N{t)/Nma.) 
is the density-dependent growth rate of all bacteria, and 
Nmax is the holding capacity, i.e. the concentration of 
bacteria supported by the nutrient broth environment. 
We directly measured the values of five parameters, a, S, 
A, X snd Nmax, and estimate the others. The dimensions 
and relevant ranges of all parameters are given in Table 

ID 

We compared the numerical simulation results with our 
experimental results for invasion and protection of the 
lysogens (Jih::(f)) carrying the phage in Fig. |3| The re- 
sults validate our choice of theoretical model and param- 
eters. We find that the experimentally determined and 
estimated values of parameters (x = 50, a = 5 /a = 0.1, 
(5 = A/a = 0.15, and the choice of P = 0.98 and 
7 = SB{0)K/a = 0.01 [7 = 0.002]) lead to a good agree- 
ment with the experimental scenarios in Fig.|3|[a) l^lb)], 
respectively. Minor discrepancies between experiments 
and simulations are noticeable in the time-evolution of 
the concentrations of BbGm::^ and of the total BbGm. 
These are likely due to our assumption of homogeneous 
mixing of phage particles and bacterial hosts. 

The success of our model in capturing the dynamic 
behavior of the bacteria-phage system in two different 
scenarios enables us to use it to determine the condition 
of a successful phage-mediated invasion and explore cases 
that are not addressed by experiments. We first investi- 
gated the invasion criterion, the choice of parameters in 
Table which makes the invading strain A dominant in 
number over the invaded strain B. The condition of be- 
coming the predominant organism depends only on the 
phage pathology P and on the initial ratio of the concen- 
trations of bacterial strains A and B. We find that the 
final population ratio is approximately 
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TABLE I: Parameters used for the numerical simulation 
of the phage-mediated competition in B. bronchiseptica. a 
[hours" ^] is determined from the measured doubling time (77 
minutes) of Bb bacteria in mid-log phase in Fig.Q A [hours~^] 
is determined from the observed phage incubation period (6- 
12 hours), which is the time-interval between initial contact 
of the phage particles with susceptible bacteria and bacterial 
lysis. X is measured from the difference in the phage concen- 
trations between and 12 hours of (j) and Bb co-culture. 5 
[hours~^] is set to be 0.054 in our simulations. Note that we 
verified in Supporting Information that the invasion criterion 
in Eq. Q remains valid for any Q < 5 / a < 0.5. The two 
undetermined parameters P and k [[hours-CFU/ml]~^] are 
estimated by comparing the experimental results with those 
of the theoretical model and by minimizing discrepancies. 
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(Free) growth rate 


0.54 


measured 
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Spontaneous lysis rate 


< 5 < a 


measured 
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(^-induced lysis rate 


0.08 - 0.17 


measured 
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Burst size 


10 - 50 


measured 


P 


Phage pathology 


< P < 1 


estimated 


K 


Contact rate 


K > 


estimated 


^ max 


Holding capacity 


~10» 


measured 



rAB(oo) > 1) is 



rABiO) >1-P. 
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The derivation of Eq. ^ was provided in Appendix in the 
hniit 7 > 1 {kS siO) > a) and /3 > 1 (A > a), that is, 
when susceptible bacteria are rapidly infected by phage 
particles and lysed immediately after infection. We also 
verified by extensive numerical simulations the validity 
of Eq. Q in a much broader parameter range, especially 
for small 7 and /3. 

The theoretical model predicts that different initial 
phage concentrations or different contact rates have no 
effect on the steady state outcome of the phage-mediated 
competition while either can modify the kinetics of the in- 
teractions (see the inset of Fig. 0}. To validate its predic- 
tion of the model, we performed time course experiments 
with three different phage concentrations (see Fig. 0J. 
The trajectories of the total BbGm population depend 
sensitively on the initial phage concentrations during the 
intermediate time-steps (between and 24 hours) , show- 
ing a larger fall of the total BbGm for higher phage con- 
centration between 8 and 12 hours of co-culture. However 
the trajectories of Bb::0 and the total BbGm converge to 
the same steady states within statistical errors after 48 
hours regardless of initial phage concentration, indicat- 
ing that the final outcome of phage-mediated competi- 
tion does not depend on the initial phage concentration. 
This verifies the independence of the invasion criterion in 
Eq. (O on the initial phage concentration and provides 
justification for the addition of exogenous phage to the 
system in Figs. Ola) and|3Ib). 

We investigated the effect of the phage pathology on 



the amount of the phage-mediated competition. We ex- 
perimentally manipulated the phage pathology using a 
lytic phage {(pAcI) with an in-frame deletion in the cl 
repressor gene required for lysogeny [l7| that always 
lyses the host bacterium (and thus has P = 1). The 
lytic phage-mediated competition experiments shown on 
Fig. Ufa) clearly corroborate the numerical simulation. 
Co-culturing strains with this (/>Ac/ resulted in an ap- 
proximately 1,000 fold advantage to the resistant strain 
(Bb: :</()) over the susceptible strain (BbGm). The small 
fraction of BbGm that survived appear to be 0-resistant 
due to infection by the wild type (j) spontaneously released 
from Bb::0. 

Finally we investigated one of several possible resis- 
tance mechanisms of bacteria against pathogenicity of 
the phage. A simple mutation can take place in the re- 
ceptor complex of the host bacteria that renders them no 
longer susceptible to phage lacking the tropism switch- 
ing mechanism. We used a mutant phage (0Aor/5) with 
in-frame deletion in the gene, orfb, encoding the reverse 
transcriptase necessary for the tropism switching mech- 
anism. This phage can only infect bacteria bearing the 
protein pertactin. In Fig.E^b) we co-cultured a bacterial 
strain with an in-frame deletion in the gene encoding per- 
tactin (Bb AprnGm) with Bb::(^Aor/5 in the presence of 
1000 PFU/ml of the mutant phage (0Aor/5) [llllllll. 
Both strains grow without any sign of phage-mediated 




40 

Time (Hours) 

FIG. 4: Independence of the steady state outcome of 
phage-mediated competition on the initial phage concentra- 
tion. Main: In vitro experiments of the time evolution of 
Bh::<f) (open symbols connected by a solid line) and the to- 
tal BbGm (filled symbols, dotted lines) with initial exogenous 
phage concentrations of 10 (filled circles), 10^ (filled squares) 
and 10^ (filled diamonds) PFU/ml. Inset: Numerical simu- 
lations of the time evolution of Bb:;0 (solid line) and the to- 
tal BbGm with three different initial phage concentrations of 
10^ (long-dashed line), 10* (dashed line) and 10^ (dot-dashed 
line) PFU/ml. The parameters are x = 50, a — 0.1, P — 0.15, 
7 = 0.02 and P = 0.98. Note that the numerical simulation 
results of the time-evolution of Bb:(j!> and the total BbGm at 
different contact rates 7 are similar to the pattern in the inset. 
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FIG. 5: Dependence of phage-mediated competition on the 
phage pathology P is observed in in vitro experiments (sym- 
bols) and numerical simulations (lines). Presented are time- 
evolutions of (a) the co-cultured Bb::<j!> (open squares, solid 
line) and total BbGm (filled circles connected, dashed line) 
in the presence of exogenous lytic phage {(pAcI), and (b) 
the co-cultured BbAprnGm (filled circles, solid line) and 
Bb::(^Aorf5 (open rectangles, dashed line) in the presence of 
mutant phage (</>Aor/5). The parameters used for the numer- 
ical simulations are the same as in Fig. inta) except 7 = 0.01 
and P = 1 in (a) and 7 = (and thus P = 0) in (b). 



competition and maintain the initial ratio for more than 
24 hours. 

DISCUSSION In most of the systems where 
pathogen-mediated competition has been observed, there 
have been confounding factors, such as different growth 
and reproductive rates of hosts, and intertwined re- 
source and phage-mediated competitions, that have lim- 
ited quantitative assessments of pathogen-mediated com- 
petition. The bacteria-phage system described here over- 
comes these difficulties to allow the accurate estima- 
tion of the parameters, in most cases directly mea- 
suring them experimentally. Understanding the im- 
pact of the pathogen pathology (the product of co- 
evolving pathogen virulence and host resistances against 
pathogens 0, M, Ull ) on pathogen- mediated competition 
requires a system in which both host defense mechanisms 
and pathogen pathology can be manipulated. Again bac- 
teria and phage are suitable systems for these purposes. 

We utilized in vitro experiments, analytical and nu- 
merical analysis to show the existence of bacteriophage- 
mediated competition between two host bacterial strains. 
In our in vitro experiments direct competition between 
host strains was minimized by using genetically identical 
strains and culture conditions that allow for unrestricted, 
exponential growth over many generations. The two bac- 
terial strains differ only in the pathology of phage infect- 
ing them; the phage pathology (P) of the strain bearing 
the prophage (Bb::(/)) is almost zero while that of the sus- 
ceptible BbGm is close to one. As we demonstrated here, 
the strain bearing the prophage has advantage in both 



invading and resisting invasion by the phage-susceptible 
bacteria. 

Our competition studies have revealed a single deter- 
mining factor for competitive advantage. First, when 
the phage was manipulated to increase its pathology in 
the susceptible strain, the amount of competitive advan- 
tage conferred to the host carrying the prophage dramati- 
cally increased. Second, when both strains (Bh::(j)Aorf5 
and BbAprnGm) are resistant to phage {4)lS.orf'b) in- 
fection, the phage didn't confer any competitive advan- 
tage to either strain. Third, we conjectured that none of 
the other details of the system contribute to the steady 
state outcome of phage-mediated competition between 
two strains, and we demonstrated experimentally the in- 
dependence of the competition outcome on the initial 
concentration of the phage. Although different initial 
phage concentrations modified the kinetics of the inter- 
actions, they did not affect the final ratio of bacterial 
strains. 

Our theoretical model captures the dynamical behav- 
ior of the bacteria-phage system and can be extended 
to predict the steady state outcome of phage-mediated 
competition under more general conditions. We demon- 
strated that when there is a quantitative difference in the 
resistance of the two bacterial strains, the success of the 
competing strains depends only on the ratio of the ini- 
tial concentrations of two strains and the relative phage 
pathology. This conclusion leads to the following pre- 
dictions for general pathogen-mediated invasion beyond 
bacteria-phage systems: a) For a pathogen to contribute 
to the ability of a host to invade an ecological niche re- 
quires that the pathogen pathology is lower in the invad- 
ing population than in the invaded population. This is 
generally true when a disease endemic to one population 
is carried to populations that are naive to the disease, b) 
In the case of differential host resistance, one can conjec- 
ture that the final ratio of the two populations and the 
success of the invasion are determined by the following 
fraction, 

Ms(oo)=r^s(0)|^^^ (4) 

where Pa[Pb\ is the pathogen pathology. The condition 
for the success of the invasion of population A to pop- 
ulation B is that r^B(oo) > 1. In Appendix we have 
verified the validity of the generalized invasion criterion 
by analytical and numerical analysis. 

The above predictions can be naturally extended to 
pathogen-mediated invasion in nature if all the details of 
individual pathogen-host interactions can be condensed 
into a single parameter describing pathogen pathology. 
There are certain limitations to extrapolating the above 
conclusions to pathogen-mediated invasion in ecology 
and humans. First, the assumption that host popula- 
tions and pathogens are well-mixed may be better suited 
in vitro than in nature. However, the theoretical models 
and experimental manipulation of phage concentrations 
suggest that the ultimate effect of pathogen-mediated 
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competition is not dependent on the rate of contact, and 
successful pathogens are, by their nature, very efficient 
at transmission (infection-causing contact). Second, the 
infection processes are discrete, stochastic and spatial in 
nature, and might not be completely described by differ- 
ential equations. Third, pathogens can be transmitted 
without killing hosts in general infection processes, and 
the birth and death events in animal and human popu- 
lations can be markedly different from bacterial growth. 
We expect that these differences will not have a major 
impact, though. 

In conclusion our studies have important implications 
in relation to the long-term advantage of bearing mul- 
tiple pathogens. An a priori view might have been that 
the heavier pathogen load might reduce fitness relative to 
a competitor. This view is directly contradicted by the 
findings of this study, in which the advantages of bearing 
a pathogen are clear and related to its relative pathol- 
ogy on the hosts. One further implication is that the 
relative advantage conferred by the pathogen to the host 
should be related to the length of time during which they 
have coexisted and co-evolved. Over time the pathogen 
will select for more resistant hosts, but will alter its own 
virulence to maintain optimal transmission and overall 
fitness. When that pathogen is introduced into a host 
population that has not been under that selection it will 
exhibit inappropriately high virulence. This effect is ob- 
served in, and potentially explains, zoonotic diseases that 
often cause more pathology in humans than in their nat- 
urally co-evolved host. Although each pathogen may 
moderate pathology to optimize reproduction in their 
new host, humans will continue to be assaulted by new 
pathogens with greater virulence. 



APPENDIX A: EXPERIMENTAL METHODS 

A temperate bacteriophage BPP-1 and a mutant, 
which lacks the reverse transcriptase necessary for the 
tropism switching mechanisms and can only infect 
B.bronchiseptica bacteria bearing the protein pertactin 
BPP-1 Aor/5, were the kind gift of Jeff Miller TElIllll. 

B. bronchiseptica strains were grown on Bordet-Gengou 
(BG) agar plate, and incubated for three days at 37°C. 
Then 2-3 colonies of each strain were inoculated in 4 ml 
of Stainer Sholte media with supplements, and grown 
at 37°C overnight to mid log phage. For co-culture ex- 
periments, two B.bronchiseptica strains were subcultured 
together into 10ml of Stainer Sholte media at the ap- 
propriate concentrations. The co-culture was incubated 
at 37°C with continuous agitation. To determine the 
concentration of each strain in the co-culture at various 
time points, lOO/il of the co-culture was serially diluted in 
PBS and spread on six Bordet-Gengou agar plates, half 
of which were treated with gentamycin. After two days of 
incubation at 37°C, the colony forming units (CFU) on 
each plate were counted to determine the concentrations 
of each strain at consecutive time points. 

The phage-sensitivity of a single BbGm colony was 
tested by using the standard protocol The lack 

of host diversity mechanisms of (/)Aor/5 was tested at 
and 48 hours postinoculation by using the standard 
phage-titering protocol on the lawns of susceptible Bb 
and BbAprnGm. 



APPENDIX B: SPONTANEOUS PHAGE 
INDUCTION 
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FIG. 6: In vitro evidence of the spontaneous release of phage. 
The strain carrying the phage (Bh::(j>, open squares) outnum- 
bers the susceptible strain (BbGm, filled circles) without ini- 
tial exogenous phage. 



In vitro evidence of the spontaneous release of phage is 
provided in Fig.|3 The strain (Bb::(/)) carrying the phage 
and the susceptible strain (BbGm) are co-cultured with- 
out exogenous phage. The initial ratio of the strain Bb:0 
to the strain BbGm is reversed around 8-10 hours, which 
is mediated by the phage spontaneously released from the 
strain Tih:(j)- Based on this result, we use spontaneous ly- 
sis rate 6 = 0.054 for the numerical simulations. 



APPENDIX C: DERIVATION OF THE INVASION 
CRITERIA IN EQS. I© AND (Hi 

Our primary model of phage-mediated competition 
depicted in Fig. 2 is limited to the case where one 
host is perfectly phage-resistant and the other is phage- 
susceptible. However in general cases both invading and 
resident hosts can be susceptible to phage infection but 
with differential susceptibilities. Here we model the in- 
vasion of a host A endogenously and exogenously car- 
rying the phage to another host B. The hosts are char- 
acterized by the differential susceptibilities ka and kb 
against the phage, and the phage pathology Pa and Pb- 
We rescale and non-dimensionalize the variables, ij = 



7 



/,/5b(0), = S,/SBiO), I, = L,/Sb{0), = ^/Sb{0), 
nraax = ^maa;/<S's (0), T = flf, Q; = 5/ a, /? = A/o, and 
7j ~ i^jSB{Q) / a where j = A, B. Then we obtain 



ds- 

— i- = (1 - n/n„iax - lj4')Sj 

dT 
di 

-j^ = {1 - Pj)-fj<j)Sj + {1 - n/n„iax - ct)ij 

^ = x(aXl*J' "ZI'^J'^^J 



(CI) 



where n = + Sj + Ij) and j =A,B. The initial 

conditions for Eq. ^fCTl are i_B(0) = /^(O) = Ib{0) = 0, 
sb(0) = 1, sa(0) > 0, iA{0) > and 0(0) > 0. The above 
general model and Eq. HC1|I are reduced to the primary 
model depicted in Fig. 2 and Eq. (1) when Pj^ — 0, — 
and Sa{0) = 0. 

CASE I: If 0(0) = and a = 0, then the 7- 
dimensional ODE system reduces to 

= (1 - (XI *J + T'A)/n„iax)sj, 
3 

= (1 - (y^ Sj + iA)lnm.ax)'iA- 



where j — A,B and ^^(t) = Ia{t) = Ib{t) — for 
T > 0. All populations sj and will grow with the 
same growth rate and the initial ratio ^^(0) : Syi(O) : 
s_b(0) remains unchanged to be ^^(t) : sa(7") : sb^t) for 
all T > O. In other words, there will be no pathogen- 
mediated competition. 

CASE II: When 0(0) > O, we can derive the inva- 
sion criteria in Eqs. ^ and in the limit 7j — > oo and 

CASE II-A: r = Limit. 

An appropriate timescale near r = is u = r/e where 
e = The effect of the transformation cr ~ t/c is 

to magnify the neighborhood of r = 0, i.e., for a fixed 
< T ^ 1, we have cr ^ 1 as e — > 0. With the trans- 
formations (T = r/e, Sj(T; e) — Sj{a; e), ij{T; e) — ij{cr; e), 
Ijir; e) = lj{a; e), 0j(r; e) = 0j(cr; e), = jj/(3, Eq. (jnU 
become 



da 



e(l - n/nmax)sj - ij4>Sj 



(C2) 



a)ij 



dlj_ 
da 



Pj4>sj£.j - ij 



3 3 



In a regular perturbation theory [23| the solutions 
are expanded in order of e, Sj (cr;e) = X]n=o 



0((t; e) = X)n=o (''')• set e = to get 0(1) 

system. 



da 
dijfi 

da 
dljo 

da 



da 



7 = x^hfi-^£.](iioSj,Q 



(C3) 
(C4) 
(C5) 
(C6) 



with the initial conditions «B.o(0)=^yi.o(0)=^B.o(0) = 0, 

SB,o(0) = 1, SA,o(0) > 0, iA.o(O) > and 0o(O) > 0. 
By integrating Eqs. ljC3p and (|C4|) . we obtain 



■Sj,o(cr) = Sj,a{Q)Exp{- I ^j(l3oix)dx), (C7) 
Jo 



+ (l-P,)s,-o(0) / F,iy)dy 



(C8) 



where f,(7/) = £_jMy)EM- lo i3Mx)dx). Eqs. (EH) 
and IIC6II can be rewritten 



dljfi 
da 



P,s,,^{0)F,{a)~l,,^{a) 



(C9) 



Lemma 1. 0o(o') *s strictly positive for a > if 4>q{0) > 
anrf x-fj > 1- 

Proof. Let Z(r7) - 0o(^t) + xEj ''j,o(ct)- ^(0) > be- 
cause 0o(O) > 0. Because ijfl{a), Sj^Q{a), ljfl{a) and 
00 (cr) are non-negative for ct > 0, ^ = S^Cx^j ~ 
l)sj{0)Fj{a) > ff x^i > 1- Therefore Z{a) is strictly 
positive and non-decreasing for all cr > 0. Suppose now 
that there exists CTq > such that 0o(cr) = for cr > ctq. 
Then both Fj{a) and lj,oi<^) will become zero for a > ao, 
resulting in Z{a) = for cr > CTq. This contradicts 
that Z{a) is strictly positive for all cr > 0. Therefore 
00 (cr) > for all cr > 0. □ 

Lemma 2. Fj{y) is strictly positive for y > and Fj{y) 
asymptotically approaches zero as y —>■ 00. 

Proof. Strict positiveness of Fj{y) for y > follows from 
lemma 1. For the second part, we divide the integration 
in the exponent into two parts. 



y ^ ry-w ^ rV ^ 

dxS,j 4>oix) = / dxS,j^o{x) + / dxS,j(j)o{x) 

Jo Jy-w 
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where y ^ 1 and w G (0,y) must be such that (l>o{x) 
is either non-decreasing or non-increasing in the inter- 
val X G [y — w,y]. Then there exists A e [0,1] such 
t^^at ly-w dxCjMx) H^^jMy) + (1 - A)Cj0o(z/ - w)]w. 



By defining 



min^>o^o{x), dx(_j(f)o{x) > 



[y - w)^j(j)rmn and Jy_^dx^j(t)Q{x) > [X^j(t)o{y) + (] 
X)^j(/)min]w. Then we can obtain 



Fjiy) = ^j4>oiy)Exp{~- / dx^j4)Q{x)) 

Jo 

< ^j$o{y)e~^'^^^^°^y'>e-^y~^"''>^^^'"^'^ 

< J_g-(y-Au,)Cj0„„„-l 

~ Xw 

where in the third hue we used xe~^ < e^^for all x > 0. 
Asy^oo, F, (y)^0. □ 

Lemma 3. Let Gj{a) = Fj{y)dy. Gj{a) asymptoti- 
cally approaches 1 as a ^ oo. 



Proof. 



dy^]$o{y)Exp{~ / dx^j$o{x)) 

JQ 

dyH'jy)e-"^^y^ 



GA<y) 



where Hj{y) ~ dx(^j(p{x) and Hj{0) ~ 0. Using 
Hj{a) > a^j4>„un, e--f^^('") < e^'^«^'^'""' for cr > 0. As 
and Gj{a) 1. 



□ 



Using the above lemmas, we know that both Sj^o(o') 
and ljfi{a) approaches zero as u ^ cxd while keeping < 
r <C 1. In the limit of cr — > oo we obtain, using Eq. I|(y8|) 
and initial conditions, sb,o(0) — 1, iB,o(0) = ^A,o(0) = 
?B.o = 0, and 



«B,o(o') + ss,o(ct) + Zb,o(o-) 
(l-P^)sA.o(0)+^A,o(0) 



(Cll) 



where rAsiO) = iA,o(0) + sa,o(0). When sa,o(0) » 
*A,o(0), Eq. 10} is recovered, in the limit of 7^ — > cxd, 
/3 — > 00 and cr ^ cxd while keeping < t ^ 1, 



(C12) 



Moreover when Pa = 0, 7a = and sa,o(0) = 0, Eq. (0) 
is recovered, in the limit of jb ^ 00, (3 ^ 00 and cr ^ 00 
while keeping < r ^ 1, 



(C13) 



rAB{(T)^rABiO)/il-PB) 



In case II-B we will prove that these ratios in Eqs. (jC12|) 
and IjClSp remain unchanged in the limit of r = 00. 



CASE II-B: T = oo limit. 

To study the long time hmit, we go back to Eq. HC1|I . 
Using the quasi-steady state approximation of the third 
equation in Supporting Eq. (1) = = Pj(f>Sjjj/(3 — 
Ij, in the limit 7^ 00 and P —>■ 00, and xPj > 1, we 
obtain ^ = xa h + ^ ixPj " 1)7, 5, > for r » 1 
(Equality holds when a — Q and Sj(r) = 0). Because 
4>{t) is strictly positive and non-decreasing for r 1 

and ^7j > 1, = (1 - n/umax - 4''lj)sj < for r > 1. 
Therefore using Ijir) = Pj(j)jjSj{T)/ Sj^r) = Ijir) — 
for T ^ 1. In the limit of /3 cx), jj — > 00 and r 3> 1, 
Eq. HC1|I reduce to the effective three dimensional ODE 



^ = (1 

dr n. 



d<p \ ^ 

d7 = ^"^^ 



Note that the first equation for ij{T) is independent of 
0(r). Now if a > 1, ^ < for T > 1 and i.j{oo) = 0, 
which means that both A and B populations go extinct. 
Otherwise if < < 1, both iA{T) and ib(t) will grow 
with the same growth rate. Thus the ratio rA_B(cr), de- 
termined in the limit of cr —> 00 while keeping < r ^ 1 , 
remains unchanged in the limit r ^ 1. 

Case III. If </)(0) = and a > 0, this is equivalent to 
case II. Suppose that 4>{t) > when r > Tmm where Tmin 
is the earliest time when the first phage are spontaneously 
induced. By defining a new time frame t' — t — Tmin 
and rescaling all the concentrations by SB{Tmin)i case 
III becomes equivalent to case II. 



APPENDIX D: NUMERICAL INVESTIGATION 
OF THE INVASION CRITERIA FROM EQS. ^ 
AND dH) 

The invasion criteria in Eqs. (O and are exact in the 
limit of large infection-induced lysis rate /3 and contact 
rate 7 with restrictions on xPj > 1 and on the spon- 
taneous lysis rate < a < 1. In order to investigate 
their validity for small f3 and 7, we performed numeri- 
cal simulations. First, the linear relationship in Eq. ^ 
between the phage pathology P and rA_B(0)/rA_B(oo) is 
validated by extensive numerical calculations with 2000 
parameter sets where all parameters are selected uni- 
formly from the biologically relevant intervals (see Fig. [7| 
for detailed information). Note that xP > 1 is used 
for numerical calculations. When 7 and (3 are relatively 
large, i.e., 0.1 < 7,/? < 10, all data points fall into 
the linear line r ab{0)/ tab {00) = 1 — P as illustrated 
in Fig. [7| When < 7,/? < 0.1, the deviation from the 
linear relationship increases for small phage pathology P. 
Thus we conclude that the linear relationship in Eq. (2), 
TAB(fi)/rABioo) = 1 — P, is robust to parameter varia- 
tions and valid for small 7 and f3. 
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Phage pathology P 

FIG. 7: Numerical verification of the invasion criterion in 
Eq. 0. A thick solid line is the prediction from Eq. ||5J. 
Tab (0) /r AS {oo) was numerically evaluated by solving Eq. 
with 2000 sets of parameters chosen uniformly in the inter- 
vals < P < 1 for phage pathology, 1/ P < x < 100 for burst 
size, < a < 0.5 for normalized spontaneous induction rate, 
< Ia{Q),4>{Q) < 105*3(0) for the initial concentrations of 
infected bacteria A and phage with respect to the initial con- 
centration of susceptible bacteria B. Filled circles represent 
the data from 1000 sets of parameters with relatively large 7 
and /3 (0.1 < 7,/3 < 10). Open circles are from another 1000 
sets of parameters with small 7 and (3 (0 < 7,/3 < 0.1). 



Second, we also validate the generalized invasion crite- 
rion from Eq. numerically with diverse sets of param- 
eters. Fig.|Slshows that the linear relationship in Eq. Q 
between r^s(0)/r^s(cx3) and (1 — Pa) /{I — Pb) is ro- 
bust against parameter variations. Note that we impose 
restrictions on xPj > 1 ^^id 3^(0) S> ia{^) in the nu- 



merical calculations. However the linear relationship in 
Eq. becomes inaccurate when the pathogen is more 
virulent on the invading population A than on the resi- 
dent population B, i.e., when Pa is large and Pb is small. 




(i-PbV(i-Pa) 

FIG. 8: Numerical verification of the generalized invasion 
criterion in Eq. I^J. rAB{Q)/rAB{oo) was numerically evalu- 
ated by solving Eq. ICH with 2000 sets of parameters chosen 
uniformly in the intervals < Pa , Pb < 1 for phage patholo- 
gies on the host A and B, l/min{PA,PB} < X < 100 for 
burst size, < a < 0.5 for normalized spontaneous induction 
rate, IO-^Sb(O) < Sa(0) < WSb{0) and < 7a(0), 0(0) < 
^"■^^^(O) for the initial concentrations of susceptible and in- 
fected bacteria A and phage. Filled circles represent the data 
from 1000 sets of parameters with relatively large 7^ and /3 
(0.1 < 7j, /3 < 10). Open circles are from another 1000 sets 
of parameters with small 7^ and /3 (0 < 7j, /3 < 0.1). 
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